Direct photons in Au+Au collisions measured with 
^ : the PHENIX detector at RHIC 
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Abstract. A major goal of experiments in heavy-ion physics is the characterization of the 
O ' quark gluon plasma (QGP) produced in the collision of heavy ions at high energy. Direct photons 

3 , are a particularly good probe of the produced medium because they do not interact strongly 

and so can escape the medium unmodified, carrying information about when the photon was 

produced. It is expected that direct photon contributions from different sources (QGP radiation, 

hard scattering, hadron gas radiation) dominate at different transverse momentum ranges. Low 

^ ' momentum direct photons are dominated by thermal radiation (both from the QGP and hadron 

0^ , gas), while high momentum direct photons dominantly come from hard parton scatterings in 

r*~ ■ the initial collision. We present a summary of techniques to measure direct photons with 

v^ ' the PHENIX detector, with a focus on low momentum direct photons through their external 

conversion to dilepton pairs. 



1. Introduction 

Experiments at the Relativistic Heavy Ion CoUider (RHIC) at BNL have estabhshed evidence 
that a quark-gluon plasma has been produced in heavy-ion collisions at ^/sJm = 200 GeV [1]. 
5^ I One signature amongst many is that the system appears to behave like a liquid with partonic 

degrees of freedom, as evidenced by elliptic flow measurements and number of constituent quark 
scaling [2]. Another observation was the suppression of high momentum probes, indicating that 
the medium is not completely transparant to colored objects [3]. Recent PHENIX measurements 
of the temperature of the plamsa also indicate a temperature well above the expected critical 
temperature needed to form a QGP predicted by lattice QCD calculations [4]. A major goal of 
the RHIC program is to now study and characterize this extreme state of matter. 

In this paper we focus on photons as probes of the QGP. Photons do not interact strongly with 
the medium and so escape virtually unmodfied, carrying information about each stage of the 
collisions where the photon was created. We are interested in direct photons, which are defined 
as all photons that do not originate from hadronic decays. This includes thermal and prompt 
photons, along with bremsstrahlung and jet-conversion photons. We mention two important 
pieces of information (although there are many others) that can be extracted from the study of 
direct photons. One is the temperature of the plasma, which can be accessed by the shape of 
the photon distribution. Another is the elliptic flow of direct photons. Theory calculations [5] 
indicate that the elliptic flow of thermal photons is quite sensitive to the thermalization time, 
To, of the plasma. It is expected that a small thermalization time leads to a small V2, as the 
proportion of photons coming from the QGP phase is larger for earlier thermalization times. 



The shape of the V2 as a function of pT also changes drasticaUy depending on the assumed tq. 
Thermal photons are expected to dominate at low momentum, and so we focus on measuring 
direct photons below 5 GeV. 

The measurement of low momentum direct photons is notoriously difficult due to very large 
background from hadron (mostly vr*^) decays. This fact, combined with the limited energy 
resolution of the electromagnetic calorimeter (EMCal) at low energy, requires an alternate 
method to simply measuring photons that directly deposit their energy in the EMCal. Thus 
it is advantageous to measure dilepton pairs originating from the conversion of photons. The 
momentum resolution in the tracking systems for charged particles improves as px decreases. 
We consider two sources of dilepton pairs from photons, internal and external conversions. 
Internal conversions occur from processes that happen to produce a virtual photon, rather 
than a real photon, which will decay into a low mass e^e^ pair [4]. External conversions 
occur when a photon interacts with material in the experimental apperature and converts into 
the dilepton pair. The external conversion analysis is complementary to the virtual (internal) 
conversion analysis and is a good cross-check of the validy of the procedure and is one of the 
main motivations for undertaking the external conversion analysis. 

2. Experimental Apparatus 

A detailed description of the PHENIX detector subsystems can be found in the reference [6]. 
Photons are measured in the central arms of the PHENIX spectrometer, each of which cover it/2 
in azimuth and total cover a psuedorapidity range of r/ < |0.35|. Charged particles are tracked 
in layers of multi-wire proportional chambers, calculating the bend in the magnetic field and 
hence the momentum of the track [7]. A combination of the Ring-Imaging Cherenkov Dector 
(RICH) [8] and the electromagnetic calorimeter (EMCal) [9] assist in electron identification. 
There are forward and backward beam-beam counters (BBCs) at 3.1 < \rj\ < 3.9, utilized as 
global detectors to determine both centrality and the vertex of the collision along the beamline 
(the z direction). A dedicated reaction plane detector consisting of plastic scintillator paddles 
was also installed for the 2007 run [10], which reside at 1.0 < \rj\ < 2.8. 

3. Analysis Techniques 

The details of the internal conversion analysis can be found in the reference [4]. The basic 
idea of this technique is that low mass e^e^ pairs are formed when a virtual direct photon is 
produced in the system. Many of these pairings exist from pure combinatorics and have no 
physical correlation. These are removed by studying pairs formed in mixed events. Further 
sources of correlated background pairs must also be removed (such as external conversions, 
cross pairs, jet pairs). These can be removed by studying like-sign foreground distributions. 
Once all background is removed, the yield of pairs is compared to a cocktail of expected 
hadronic sources. An excess of pairs is observed in Au+Au collisions and is interpreted as 
originating from virtual photons. The yield of real photons can be calculated from the pair rate 
observed. Upon calculating the cross-section of direct photons, an excess above the expected 
rate is seen below pr < 3 GeV. This excess is exponential in shape and is fit to extract a 
temperature parameter. The temperature measured in the paper for central Au+Au collisions 
is T = 221 lb 19{stat) it 19{sys) MeV, which is above the expected critical temperature needed 
for the phase transition. This is an important result and it is worthwhile to seek independent 
measurements of the same quantity as a cross check. The external conversion channel provides 
us with this cross check. 

The details of the external conversion analysis are described here. We start by identifying 
single electrons, which will be paired later. The electron id detectors used are the RICH and 
EMCal using common PHENIX electron id cuts. It is required that the track is associated with 



> 2 phototubes in the RICH. It is also required that the energy to momentum ratio, E/p > 0.5. 
The reconstructed momentum of the track must be 0.2 < pj- < 20 GeV. 

We need to identify pairs that actually come from an external photon conversion. During 
the 2007 RHIC run, the Hadron Blind Detector (HBD) [11] was installed for its engineering 
run. The backplane of the detector is fairly thick (about 4%Xo) and so is a significant localized 
source of conversions. Hence we choose to focus on conversions occurring in this backplane. It 
is located at a radius of about 60 cm from the interaction point. This significant conversion 
radius is the key to indentifying these photon conversions. 

PHENIX does all of its tracking outside the magnetic field. This requires an assumption 
about the origin of the track to determine the bend in the field and calculate the momentum 
of the particle. A natural assumption is that tracks originate from the event vertex. But if we 
consider electrons coming from a radius of 60 cm, then this assumption is wrong and leads to 
a misreconstruction of the electron track. This misreconstruction leads to an artificial opening 
angle of the pair, which leads to an apparent mass of the pair. This apparent mass is proportional 
to the radius at which the conversion occurs (the controlling parameter is how much field the 
particle missed, even though it was assumed to pass through the entire field), causing a slight 
over-estimate of the momentum. The apparent mass can therefore be used to help identify 
conversions at the HBD backplane. Simualtion studies have shown that we expect an invariant 
mass peak of the dilepton pairs to be at 12 MeV with the Run 7 magnetic field configuration, 
see the open symbols in Fig. 1. This same peak can be seen in the data. Fig. 2. 
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Figure 1. The e~^e~ pair invariant 

mass distriubtion for photon conversions at 
different radii is plotted from GEANT based 
Monte Carlo simulations. The source of the 
conversion is represented by various shapes as 
shown in the legend. Open symbols represent 
the normally reconstructed mass, with the 
closed symbols as the ATM calculated mass. 




Figure 2. The application of the pair cuts 
on the raw pair invariant mass distribution 
shown for one particular px bin in real data. 
The left side shows the normal track model 
and the right side the ATM. Each color 
represents the application of more cuts as 
shown in the legend. 



Of course it is desirable to have the correct momentum of the converted photon, not a 
misreconstructed momentum. Therefore we correct for this mis-measurement with an alternate 
track model (ATM) assumption. Rather than assuming that tracks originate from the event 
vertex, it is assumed that all the tracks originate from a radius of 60 cm. The momentum of the 
reconstructed track (under the normal track model assumption) is reparameterized for the ATM 
assumption. This is studied through Monte Carlo simulations. The conversion pairs are then 
reconstructed properly under this ATM assumption, and the invariant mass moves from peaking 
at 12 MeV towards having a peak at MeV. All other particles originating from the event vertex 



will now be misreconstructed, but this is a desired effect. Pairs from particles originating from 
the event vertex will be shifted up in mass. This is especially important for Dalitz decays of 
vr'^'s, since the HBD conversion peak sits on the side of the Dalitz peak (see the left hand side of 
Fig. 2, where the first peak seen in the black historam is the Dalitz peak and the second peak 
from conversions at the HBD). We cut on how the mass moves under the two track models, 
allowing us to indentify converted photons. It is required that the ATM calculated mass is less 
than 4.5 MeV and the normally calculated mass is between 10 < mass < 15 MeV. 

In addition to these mass cuts, we also cut on the separation of the tracks at the HBD 
backplane calculated with the corrected ATM momenta. Conversion pairs should have no 
opening angle at the point of conversion. Therefore we require that the tracks are close in (j) 
and z. These cuts allow us to get a very clean photon id at low momentum. The effect of these 
pair cuts on the raw e'^e'^ pair invariant mass distribution is shown in Fig. 2. An unfortunate 
detector problem offers a good opportunity for an internal cross-check of the analysis, namely the 
removal of half of the HBD for repair during the run. And so we can check the misidentification 
rate of the photons by analyzing the pairs we id in the arms serparately and comparing the 
number of photons that we find to convert in the arms. It is observed that the misidentification 
of converted photons is less than 3%. The combinatorial background is insignificant in this very 
low invarant mass range (mass less than 20 MeV). 

The final goal of the measurement is the direct photon V2 (at the time of writing this analysis 
is still a work in progress and so only an analysis strategy and current status is presented). 
We need three basic pieces of information to derive the direct photon V2- We measure the V2 
of all photons (we cannot a priori distinguish the origin of a photon) and statistically subtract 
background sources to the V2- Therefore, we measure the inclusive photon V2- In addition 
to the inclusive photon V2, we need to also know the vr'^ decay photon V2 (since most of the 
inclusive photons come from vr'^ decays). This is the only hadronic V2 source that we measure. 
Contributions to the V2 from other (less significant) sources of hadrons will be calculated from 
KEt scaling, where KEt = rriT — rriQ. We also need to know the relative fraction of direct 
photons compared to hadronic sources or R^ = jv*" /A^ ", the ratio of inclusive photons to 
hadronic decay photons. Then the direct photon V2 can be calculated from Eqn. 1. 

]\Tincl„,inc ]\tBG„,BG n „,ind n,BG 

dir ^ iV V2 - ^V V2 ^ n^V2 " V2 ,-.-. 
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The inclusive photon V2 is measured by taking the identified external photon conversions using 
the aforementioned method and correlating the angle of emission, (/>, relative to the reaction plane 
angle, ^, in azimuth. This (p — '^ distribution is binned in px- Correlation functions (C.F.) are 
constructed from these distributions. We perform a Fourier decomposition of the correlation 
functions, following Eqn. 2. The correlation functions are fit with Eqn. 2 and the elliptic flow 
strength, V2, is extracted in each px bin. Higher order Fourier terms in Eqn. 2 are assumed 
to be negligible and are not included. Sine terms are zero by symmetry. We also calculate V2 
as {cos{2(4> — ^))), the mean projection method. This yields consistent results with the fit and 
the average of the two is taken for the result presented here. The extracted V2 is corrected 
by the reaction plane resolution, determined by comparing correlations between the north and 
south detectors. The major source of systematic error on the points is from the reaction plane 
determination, which is roughly 10%. 

C.F. = B[l + 2v2Cos{2{(/)-'^))] (2) 

The ratio R^ is measured using the double ratio shown in Eqn. 3. It is defined as the ratio of 
the yield of inclusive photons to the yield of hadronic decay photons. A ratio above one indicates 
a signal of direct photons. The terms that go into the numerator of the double ratio, along with 



the correction factors, are shown in Eqns. 4 and 5. Note that ah the terms are measured as a 
function of the converted photon pT- The numerator of Eqn. 3 is purely from data and consists 
of the ratio of the inclusive photon yield (the converted photons) to the yield of those photons 
that we tag as coming from vr*^ decays. It is necessary to measure the tt^ contribution since 
that is the largest source of hadronic decay photons. This is done by reconstructing tt'^'s by 
combining the converted photons with photons we measure in the EMCal. To identify a photon 
in the EMCal we cut on the shower shape, comparing the shower shape of the photon candidate 
with the expected shower shape for a photon, and require a minimum energy deposition of 500 
MeV. The combinatorial background is removed with a mixed event technique and after this 
subtraction, we integrate the mass distrubution around 2a of the vr'^ peak. In this way we 
(statistically) tag some of our inclusive photons as coming from a vr'' decay. 






N'^'^ipT) = Cepa^rapa^n''"'\pT) (4) 

Nf'^'ipT) = Ceparrapaire^fY\pT) (5) 

In the Eqns. 3, 4, and 5, 7 represents the yield of photons that nature produces, N represents 
the value we actually measure. The correction factors for the inclusive photon yield are shown 
in Eqn. 4. epair and apair are the e^e~ pair efficiency and acceptance in the PHENIX detector 
respectively. There is also some factor corresponding to the probability of the photon externally 
converting, denoted by c. These corrections apply to both the inclusive and the 7r° tagged 
photon samples. For the vr*^ tagged sample we additionaly have to correct for the reconstruction 
efficiency, e-y, and apply an acceptance correction, f, for reconstructing the unconverted decay 
photon in the EMCal. The f factor is a conditional acceptance correction and is defined as 
the probability of getting the unconverted photon in the EMCal, given that we already have 
the e"^e^ pair from the other photon. The e^ accounts for occupancy effects in the detector 
by embedding simulated tt^ decays into real events. These corrections are in progress and are 
being studied with full GEANT3 based Monte Carlo simulations. We account for sources of 
decay photons other than the ir^ in the denominator of Eqn. 3 determined from simulations. 
Main sources of decay photons are decays of mesons such as tt'^, 77, and ??', and are included in a 
photon cocktail. The photon distributions from these decays are calculated in a decay generator, 
which realistically handles the kinematics of the decay. The measured tt^ pt distribution from 
published data is fit with a Hagedorn function as is used as input to the photon cocktail for 
the TT*^ Pt shape. The shape of the pT distributions for the other hadrons is determined by rriT 
scaling, where the pT is replaced by rriT = yjfn + pt, with m the mass of the hadron. The 
relative hadron/vr*^ ratios are also taken from actual measurement and input into the cocktail. 
Then the denominator of the double ratio can be calculated, which is the ratio of the yield of 
all hadronic decay photons to the yield of decay photons coming from n^s. 

The vr'^ tagging technique will also be extended, allowing us to calculate the ir^ V2 in a manner 
consistent with the rest of the analysis. 

4. Results and Conclusions 

This analysis is still a work in progress at the time of writing, therefore only the methodology 
and current status have been discussed. A preliminary result has been obtained for the inclusive 
photon W2, shown in Fig 3, for three centrality bins in a range 0.4 < px < 3.5 GeV. From left 



to right in the figure are the — 20%, 20 — 40%, and 40 — 60% centrahty bins. The major 
source of systematic error on the points comes from the reaction plane determination (about 
10%), with photon id and calculation method contributing only a few percent. The systematic 
error associated with the reaction plane resolution is common to both the inclusive photon V2 
and the tt^ V2, thus when the final systematic error is propagated to the full direct photon V2 
result, the reaction plane systematic will only come in once. These results are compared to 
published data [12] in which photons are directly measured in the EMCal. The agreement is 
good, indicating the validity of this method. 
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Figure 3. The preliminary measured inclu- 
sive photon V2 via the external conversion 
method. This is shown in three centrality 
bins. From left to right the bins are — 20%, 
20 - 40%, and 40 - 60% centrality. 
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Figure 4. The same points as in 3, but with 
published data overlayed in the open black 
dries with the red systematic error boxes [12] 



In summary, we have presented the current status and the analysis method for measuring 
the elliptic flow of direct photons via their external conversion to dilepton pairs in the detector 
material. A comparison of the measured inclusive photon V2 with an earlier PHENIX result 
obtained using the EMCal shows the reliability of this method. In a next step, we will subtract 
the V2 of decay photons and calculate the direct photon V2, using the external conversion method. 
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